Particle multiplicities and particle ratios in excluded volume model 
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One of the most surprising results is to find that a consistent description of all the experi- 
mental results on particle multiplicities and particle ratios obtained from the lowest AGS to 
the highest RHIC energies is possible within the framework of a thermal statistical model. 
We propose here a thermodynamically consistent excluded- volume model involving an inter- 
acting multi-component hadron gas. We find that the energy dependence of the total mul- 
tiplicities of strange and non-strange hadrons obtained in this model agrees closely with the 
experimental results. It indicates that the freeze out volume of the fireball is uniformly the 
same for all the particles. We have also compared the variation of the particle ratios such as 
(K+)/{*+), {K~) /(*-), K-/K+, p/p, A/A, B/S, fi/fi, (A)/<7r), (S-)/<7r), (fi + fi)/<7r) 
and ($)/(7r) with respect to the center-of-mass energy as predicted by our model with the 
recent experimental data. 
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I. INTRODUCTION 



In ultra-relativistic heavy-ion collision, we expect that a matter with high energy density will 
be produced and finally a colour-deconfined matter is likely to be produced. Such a matter, better 
known as quark-gluon plasma (QGP), is only a transient state of the system because evolution back 
to ordinary hadronic matter begins immediately [H, 0, 0). Consequently we face a difficult question 
how to diagnose and probe the plasma state. One of the most potential probe for detecting QGP 
formation is the observation of the abundance of the strange particles. This is based on the proposal 
that in the baryon rich plasma, the strange quark-antiquark pairs ss would be more abundant than 
u, d quarks because of the existence of Pauli blocking in light quark creation. The threshold for ss 
pair production (~ 300 MeV) in QGP is much lower than that for KK production (~ 980 MeV) in 
HG. Based on these arguments, Rafelski and Miiller have suggested that a large 

strange quark present in QGP would automatically yield an enhanced particle production in the HG 
resulting after hadronization of QGP. However, one important assumption in the theory is that the 
HG formed after QGP phase does not find sufficient time to achieve chemical equilibrium, otherwise 
the equilibrated HG will not retain the memory regarding the primordial phase that might have 
existed earlier in the evolution process. However, in the analysis of strangeness enhancement, one 
should have the precise knowledge not only of the equation of state (EOS) of QGP alone but of 
HG also which gives the background contribution in this case. 

In ultra-relativistic nucleus-nucleus collisions, a very high density of matter exists over an 
extended region, and it is often called a 'fireball'. The physical variables characterizing a fireball 
are the energy density e, the baryon density ub and the volume. Thus one of the central problems 
in the study of high energy collisions lies in deducing the state of the system by determining 
temperature and baryon chemical potential [is (which in turn determine e and ns) existing in 
the fireball from the observed final particle multiplicities etc. The freeze-out stage of the system 
when the particles fly towards the detectors without further interactions is directly connected to 
the observed multiplicity distributions of various particles. Here we should emphasize that the 
fireball goes through two types of freeze-out stages. First it suffers a chemical freeze-out when 
inelastic collisions in the fireball do not occur anymore. Later when elastic collisions also stop 
in the fireball, the stage is called as thermal freeze-out. Abundances of the particles and their 
ratios provide important information regarding the chemical equilibrium occurring in the fireball 
before the thermal equilibrium. Moreover, chemical equilibrium in the hot, dense HG removes any 
memory existing in the fireball regarding a primordial phase transition. So the thermodynamical 
properties of the fireball do not throw any light on the existence of QGP before hadronization. 
We first connect the thermodynamical properties of the fireball to an appropriate EOS for the hot 
and dense hadron gas and finally we deduce the chemical freeze-out conditions of the thermal HG 
fireball formed in the ultra-relativistic heavy-ion collisions. In the last few years, one has observed 
a very surprising result. It has been shown that a consistent and appropriate description of all the 
experimental results on the particle multiplicities and particle ratios from the lowest SIS to the 
highest RHIC energies is available within the framework of a thermal statistical model as applied 
to a hot and dense HG phase. Recently we have given a model based on the geometrical excluded 
volume correction which describes suitably the thermodynamical quantities of a hot and dense 
HG [HI]. We have further used this prescription to determine the freeze-out volume of the fireball 
as well as pions and nucleon density [l2|]. We find that the densities are reproduced well by their 



HBT experimental data. 



In this paper we attempt to explain the experimental data on strange particle multiplicities, 
particle ratios in the HG scenario using our excluded volume model discussed in 11, 1^]. Since many 
131 1 have been made by various authors to explain particle multiplicities using different 



attempts 



lq . 1 171 . Il8l ] , it will be worthwhile to compare chemical freeze-out parameters 



HG models 

(i.e., T, fig) °f the fireball as determined in these models. Finally we show our curves for strange 
particle multiplicities and particle ratios and compare them with the experimental data. We also 



show the predictions of the other important models existing in the literature [la . Il9l . l2Qf ] . It is 
indeed encouraging to notice that simple thermal models can explain well the experimental data on 
particle multiplicities and the strangeness enhancement. But the purpose here is to demonstrate 
that the whole exercise crucially depends on the EOS used to determine the thermodynamic state 
of the fireball of the hot and dense HG. 



II. METHOD OF CALCULATION 



We have formulated a new thermodynamically consistent excluded-volume model HI] for the 
description of hot and dense hadron gas (HG). Here we are describing in brief for the sake of 
completeness. This model assumes that the baryon of i th species have an eigen volume V{. If 
R = n \ x Vi be the fraction of occupied volume, then the number density n\ x of i th baryon can 
be written as: 

B /? 

nf = {\-R)I i \-I i \\ — , (1) 

where Aj is the fugacity of i th baryons and Ij is the following expression containing modified Bessel 
function of second kind 

with gi is the spin-isospin degeneracy factor. Eq. (1) can be rewritten in the form 

R=(l-R)J2n^V i -J2^V i X i ^, (3) 

i i 

with n° = Ii Aj. Taking R° = ^ Xi, where X% = Ii Aj V{ and putting dR/dXi = 0, we get 

F?° 

R = R = — — — 77. (4) 

Thus we can write Eq. (3) in the form 

R = R + £IR, (5) 

with is given by 

d 



n 



Now by following Neumann iteration method Eq. (3) can be cast in the form 



R = R + nR + n 2 R + n 3 R + ■■■ (7) 

Retaining terms upto second order only, the expression for R can be written as 

Finally by calculating the values of R and dR/dXi, one can calculate the value of particle 
number density by using Eq. (1). 

We have considered here a hot and dense hadron gas with baryonic and mesonic resonances 
having masses up to 2 GeV/c 2 . In order to conserve strangeness quantum number, we have used 
the criteria of net strangeness number density equal to zero. In all the above calculations we have 
also considered that all the mesons behave as point-like particles. Furthermore, we have taken 
equal volume V = 4-7T r 3 /3 for all baryons with a hard-core radius r = 0.8 fm. 

We consider all hadronic resonances having well defined masses i.e., their decay widths are 
small. All hadronic resonances decay rapidly in strong decays after freeze-out and thus contribute 
to the stable particle abundances. Some heavy resonances may decay in cascades. This has been 
implemented in the calculation by considering all decays proceeding sequentially from the heaviest 
to lightest particles. As a consequence of this, the light particles get contributions from the heaviest 
particles. Thus it has the form : 

m = b 2 ^\ ■ ■ ■ b N b N ^i n N , 



where b^k-i combines the branching ratios for the A; — > A; — 1 decay 21| with the appropriate 
Clebsch-Gordan Coefficients. The latter accounts for the isospin symmetry in strong decays and 
allows us to treat separately the different charged states of isospin multiplets of particles. In 
order to calculate contributions from some heavy resonance decays we are forced to take some 
approximations. For example, in most of the cases there are several decay channels. In our approach 
we discard all decays with branching ratios less than 2%. In addition if the decay channels are 
classified as dominant, large seen or possibly seen, we take into account the dominant channel 
only. If two or more than two channels are described as equally important, we take all of them 
with the same weight. For example, /o(980) decays into it (according to [21] this is the dominant 



channel) and KK (according to 22] this is the seen channel). In our approach, as a rule stated 
above, we include only the process /o(980) — > ir. Similarly ao(1450) has three decay channels: t]tt 
(seen), tttj (958) (seen) and KK (again seen). In this case we include all the three decay channels 
with equal weight 0.33 (branching ratios). Of course, this procedure is not unique and may vary 
from author to author. Another problem that we face in this calculation is that in some decay 
channels the branching ratios are not given exactly but a range of values are given and the sum of 
the branching ratios may differ significantly from 1.0. In this case we take the mean values of the 
branching ratios. Since we require that their sum should be properly normalized, we are forced to 
rescale all the mean values in such a way that their sum is indeed 1.0. Since the experimental data 



TABLE I: List of particle ratios used to fit different models [l^ . 
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on the resonances involve a lot of uncertainties in decay width as well as in branching ratios, we 
adopt the above procedure in calculating the contributions of resonances towards a given particle 
multiplicity. Unfortunately, a better and established alternative for such calculation does not exist 
in the literature. 



III. RESULTS AND DISCUSSIONS 

In order to find out energy dependence of particle ratios and multiplicities we first find energy 
dependence of chemical freeze-out temperature and baryon chemical potential by fitting the different 
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TABLE II: (T, //g) values in MeV obtained by fitting the particle ratios using different models. 





Ideal 


gas model 


Rischke model 


Present model 


\/SNN 

(GeV) 
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S 2 
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S 2 


T 
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2.70 
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particle ratios from SIS to RHIC energies using ideal HG model, Rischke model as well as present 
model. In Table 1, we have given the list of particle multiplicity ratios [13| which had been 
used to calculate T, mb values at freeze-out for different energies. Table 2 summarizes the values 
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FIG. 1: (a,b,c). Variations of the 
strange meson to non-strange meson ratios 
(K+)/{ir + ), (K-)/{ir-) and anti-kaon to 
kaon ratio K~ /K + with respect to ccntcr- 
of-mass energy ^snn- Solid curve (C), 
dashed curve (D) and dash dotted curve (E) 
show the predictions of the present, ideal 
HG and Rischke model, respectively. B rep- 
resents experimental points [29| . 



of T and jug at different center-of-mass energies in different models. We find that T and hb 
obtained in our model can be parameterized in terms of center-of-mass energy by using following 




FIG. 2: (a,b,c,d). Variations of the anti-baryon to baryon ratios p/p, A/A, S/S and fi/fi with respect to 
center-of-mass energy. B represents experimental data whereas curve C, D, and E show the results from 
present model, ideal HG model and Rischke model, respectively. 





FIG. 3: (a,b,c,d). Variations of the strange baryon to non-strange meson ratios (A)/(ir), (E) / (tt) and (f2 + 
0) /(it), ($)/(7t) with the center-of-mass energy \J Snn- In this figure curves B represents experimental 
data [30( and curves C, D and E depict the predictions of present model, ideal HG model and Rischke 
model, respectively. 




FIG. 4: Variations of total multiplicities of A, (f2 + f2), E - , A, <I>, K~ , K + and 7r + with respect to center- 
of-mass energy predicted by our present model. Experimental data measured in central Au+Au/Pb+Pb 
collisions [37|-[53( have also been shown for comparison. In this figure A represents multiplicity of A (scaled 
by a factor 0.02), B that of U, + Cl (scaled by 0.2) C that of S~ (scaled by 0.1) and D that of A (scaled by 
0.02). Similarly E, F, G and H represent the multiplicities of $, K~, K + , ir + mesons, respectively. 



equations [23, U& 0, 01 



MB 



1 + b tJs nn ' 
T = c -de~ f v 7 ^. 



(9) 



Here the values of the parameters as arising from the best fit are a = 1.254 db 0.086 GeV, b = 
0.261 ± 0.034 GeV" 1 and c = 168.49 ± 4.4 MeV, d = 171.63 ± 21.75 MeV / = 0.21 ± 0.039 GeV. We 
compare our values of the above parameters with the values of the parameters a = 1.308 ± 0.028 
GeV, b = 0.273 ± 0.008 GeV" 1 and c = 172.3 ± 2.8 MeV, d = 149.5 ± 5.7 MeV / = 0.20 db 0.03 
GeV obtained by Cleymans et al., in ideal HG prescription [25 



261 ] . We find that the difference 



in the values of the parameters as obtained by us and by the earlier authors does not appear very 



significant. Recently it has become abundantly clear that the heavy-ion collision energy plays an 
important role in determining the properties of the final state hadrons. We find that the extracted 
temperature generally increases rapidly whereas baryon chemical potential monotonically decreases 
with the collision energy. In general, extracted freeze-out values of these parameters in our model 
lie close to the ideal gas values. The parameter values in Rischke model are found to differ in 
comparison to our values. 

Particle ratios and multiplicities are then calculated using these T and \ib- In order to cal- 
culate the multiplicity of hadrons, we first determine the total freeze-out volume by dividing the 
experimentally measured multiplicities of K + meson at different energies with its number density 
as calculated in the present model at different center-of-mass energies. We assume that the fireball 
after expansion, goes through the stage of chemical equilibrium among its various constituents 
and the freeze-out volume of the fireball for all types of hadrons at the time of the homogeneous 
emissions of hadrons remains same for all particles. This freeze-out volume thus extracted, has 
been used to calculate multiplicities of all other hadrons from corresponding number densities at 
different center-of-mass energies. 

In Fig. l(a,b,c), we show the center-of-mass energy dependence of the ratio of strange meson 
to non-strange meson such as (K + ) / (ir + ) , (K~) and K~/K + in our present model, ideal 
hadron gas model and Rischke model. Solid, dashed and dash dotted curves represent predictions 
of our model, ideal hadron gas model and Rischke model, respectively. The {K + ) / (ir + ) ratio 
shows a peak at around 8.0 GeV of center-of-mass energy. While the ratio (K~) / (ir~) shows 
a monotonic increase with center-of-mass energy x/sjvjv an d achieves a saturation at or around 
RHIC energy. Essentially the saturation arises as the freeze-out temperatures that become almost 
constant at or around RHIC energies. The peak in the (K + ) / (ir + ) ratio is demonstrated in all the 



thermal models [27|, [28| and this successful explanation of the experimental data is one important 
feature of the thermal models. However, the peak in these models is slightly broader than that 
observed by experiments. Our thermal model shows a peak at the same energy as measured in 
the experiment and we find a good agreement between our thermal model prediction and the 
experimental results (23 |. Other types of models like transport model etc. are not able to explain 
the experimental feature of (K + ) / {tt + ) ratio. However, all types of thermal models provide a large 
ratio for K~ /ir~ than observed in experiments. Still the main feature of the data showing a steady 
increase is also reproduced by all thermal models. The ratio K~ /K + shows much dependence 
on the center-of-mass energy. Our thermal model prediction for the ratio K~ /K + shows a close 
agreement with the predictions of other models and also matches with the experimental results. The 
feature of the experimental curve is also in a good agreement with our thermal model calculation. 

In Fig. 2(a,b,c,d), we plot the energy dependence of the anti-baryon to baryon ratios like 
p/p, A/A, S/E and Cl/Q as given by our present model, ideal hadron gas model and Rischke 
model. Surprisingly we find that our present model and the ideal hadron gas model reproduce 
the qualitative and quantitative features of the experimental data for all of these ratios. This 
occurs because the freeze-out values of T and at different energies do not differ much in these 
two models. However, the predictions from the calculations in Rischke model do not show such 
agreement with the experimental data. As we have demonstrated in the previous chapter that 
ideal HG description does not provide a proper EOS for hot and dense hadron gas. The successful 
explanation of the data shows that a proper and realistic EOS for HG is given by our model which 



gives a thermodynamically consistent description of the hot and dense HG. 

Fig. 3(a,b,c,d) show the variation of the particle ratios {A)/(n), (3)/(7r) and (0 + 
Q)/{tt), (&) / (it) with the center-of-mass energy ^snn in all the above three models. We find 
that results in our present calculation particularly of (A)/{tt) show good agreement with the ex- 
perimental data [30(] since our curve reproduces the main feature of the data. Ideal hadron gas 
model also gives a satisfactory explanation but results lie far below the experimental points. This 
again supports our earlier claim that a suitable description for HG is given by the EOS in our 
model. Main problem comes when we compare predictions from all these models for multi-strange 
particle productions e.g.,(H) / (tt) and (fi + f2)/(7r), ($}/(7r) as shown by Fig. 3(b,c,d). We find 
that in all such cases, thermal model predictions do not match with the experimental data and in 
fact theoretical curves lie far above the curves given by the experimental data. Here S is ssu (or 
ssd), il is sss and is given by ss quark combinations. These results not only signify the failure of 
all types of thermal models but also emphasize the need for going beyond such approaches for the 
cases of multi-strange hadrons. The strangeness enhancement signals QGP formation but in these 
cases, thermal models yield much larger str ang eness than observed by the experiments. However, 



it has been pointed out by several authors [311 , 132 . l33l . |34| . l35l . l36j ] that such anomalous results are 



compatible with quark-coalescence models invoking a QGP formation. 

In Fig. 4, we show the center-of-mass energy dependence of multiplicities of hadrons A, (Q + 
fj), E~, A, <£, K~ , K + and tt + as predicted by calculations in our present model. Curves A, 
B, C and D show the multiplicities of A, $7 + Cl, and A baryons scaled by factors of 0.02, 
0.2, 0.1 and 0.02, respectively. Curves E, F, G and H depict the multiplicities of K~ , K + 
and 7r + mesons, respectively. We have also shown here experimental results measured in central 



Au+ Au/Pb+Pb collisions [37|-[53j] for comparison. Here we have extracted freeze-out volume of 
the fireball from the calculated number density of K + and comparing it with total multiplicity of 
K + experimental data. We use the same volume for all other particles. We observe an excellent 
agreement between the theoretical predictions by our present thermal model and the experimental 
data for the total multiplicities of n + , K + , K~ , A, A etc. However, thermal model calculation 
slightly differs for + 0, H and <!> as compared to experimentally measured values. Thermal values 
of the multiplicities for all these particles are again larger than the experimental values. This 
analysis again suggests a new and different mechanisms for the production of these particles. One 
way out of this difficulty is to assume that these particles achieve chemical equilibrium earlier in 
the fireball when the corresponding volume is much smaller. However, this is a vexing problem 
appearing in the use of thermal models and we have to find appropriate answers to these problems. 



IV. SUMMARY AND CONCLUSIONS 



We have formulated an excluded-volume hadron gas model and used it in order to analyze the 
variations of multiplicity of various particles with respect to the center-of-mass energy ^snn from 
SIS to RHIC energies. We have also used the present model to explain the variations of some 
particle ratios again with center-of-mass energy and compared our results with the experimental 
data. A good agreement between our present model results and experimental data supports the 
claim that thermal model gives a satisfactory description of the data. 



In the past we have witnessed the success of ideal gas model in explaining the data, but we 
know that a correct description of hot, dense hadron gas can be given by a model where hard- 
core repulsive interactions are incorporated in a thermodynamically consistent way. Moreover, at 
GSI-SIS and BNL-AGS energies, the freeze-out parameters involve a much larger values of baryon 
chemical potentials and the predictions of all the excluded-volume models are quite different from 
those obtained in the ideal gas models. We have alread y sh own the utility of the present model in 
explaining various properties of hot, dense hadron gas [111 . The analysis presented here lends 
further support to our claim that the excluded-volume model obtained by us properly explain the 
multiplicities and particle ratios of various particles after chemical freeze-out. In conclusion, our 
model provides a proper and realistic EOS for a hot, dense hadron gas and it can successfully be 
used at extreme values of the temperatures and/or densities. 
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